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Abstract: Some part of the relic Dark Matter is distributed in small-scale clumps which survived struc- 
ture formation in inflation cosmological scenario. The annihilation of DM inside these clumps is a strong 
source of stable charged particles which can have a substantial density near the clump core. The stream- 
ing of the annihilation products from the clump can enhance irregularities in the galactic magnetic field. 
This can produce small scale variations in diffusion coefficient affecting propagation of Cosmic Rays. 



Introduction 

The Cosmic Ray (CR) propagation below 10 17 eV 
can be described as a resonant scattering on the 
magneto hydrodynamic turbulences (MHD) with 
the scale k equal to the particle Larmor radius r g 
in the galactic magnetic field B, k^ 1 ~ r g = 
pc/ZeB [1]. The MHD turbulences can prop- 
agate in space as Alfven waves with the veloc- 
ity v a ~ B I \/4irpH, which depends on the in- 
terstellar gas density pn and is in the order of 
10 7 cm/ s. The level of MHD turbulences is pro- 
portional to the random component of the mag- 
netic field SB/B and can be associated with fluc- 
tuations in interstellar medium (ISM) which fol- 
low a power law W(k) oc k a ~ 2 in the range of 
k = 1(T 20 - lO^cm" 1 , where a = 1/3 for the 
Kolmogorov spectrum. The waves interact with 
the CR and interstellar medium and can be en- 
hanced or damped, depending on the energy flow. 
In the self consistent approach the kinetic equation 
for the spectral density W(r, k) of MHD turbu- 
lences in spherical coordinates can be written as 
[1]: 



dW 1 d ^ ^ w dVg d 
dt r 2 dr a dr dk 



kW = (G-S)W (1) 



The G term describes the enhancement of turbu- 
lences due to streaming of CR particles and S rep- 
resents the damping. The growth of turbulences 
occurs when the CR streaming velocity v s is larger 
than the Alfven speed v a [2]. The growth depends 
on the gradient of the CR density / which can 
be obtained from the steady state solution of the 
diffusion equation (here without convection and 



reacceleration)[3]: 
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f - avriH f = -q(p,r) (2) 



, where q(r,p) is the source term, || energy losses 
and a fragmentation cross section. The diffusion 
coefficient D at resonance is related to W(k, r) as 



D(r) 



)B 2 /l2irW(k r ,r) and a high level 



of turbulences corresponds to the small coefficient 
and local confinement of particles. Since the en- 
hancement and damping strongly depend on the lo- 
cal environment, this opens a possibility for spatial 
variations in propagation parameters and therefore 
CR density. Here we consider how the dark mat- 
ter (DM) annihilation can introduce such variations 
and affect the CR propagation. 

Dark Matter annihilation in clumps 

The N-body cosmological simulations and analyt- 
ical calculations show that in the inflation scenario 
the smallest DM structures, or clumps, originate 
from primordial density fluctuations. These pri- 
mordial clumps are partially destroyed during evo- 
lution contributing to the bulk DM but 0.001-0.1 
of total relic DM can still reside in clumps, de- 
pending on initial conditions [5], The clump mass 
distribution follows n(M)dM ~ M~ 2 with the 
minimum mass M™ in - 10" 8 - 1O- 6 M de- 
fined by free streaming of DM particles after ki- 
netic decoupling [4]. The local number density 
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distribution of clumps n c i depends on the bulk 
density profile and tidal destruction in the galaxy 
[6]. The clumps density profile is probably cuspy 
Pd oc I/?- 1 - 5 - 2 - but is saturated at some criti- 
cal density p max forming a dense core r c [4]. In- 
side the clump, the DM of mass m x will annihi- 
late producing stable particles: protons, antipro- 
tons, positrons, electrons and gamma rays, which 
can be observed on top of the ordinary CR fluxes. 
The luminosity of the clump for an i-component 
is: qi(r,p) ~ {<Tv) Jj {p) p 2 max rl, where Y, is the 
yield per annihilation. For most of DM candi- 
dates the annihilation goes into fermions, predom- 
inantly quark-antiquarks pair. For example the 



despite of smaller luminosity the local density of 
produced particles can exceed the galactic average 



mSUGRA neutralino of 



100 GeV anni- 



hilating in bb will produce per annihilation at 1 
GeV : ^3 positrons or electrons, 0.3 protons or an- 
tiprotons and 8 gammas, the precise energy spec- 
trum from quarks hadronization is well measured 
in accelerator experiments. The (av) is the ther- 
mally averaged annihilation cross section which 
can be estimated from the observed relic DM den- 
sity in time of decoupling (Tdec ~ -jjf ): (<rv) w 
( 2 - 1 °'^ 3s ' 1 ), where ^ x /t 2 = 0.113±0.009 [7]. 
Nowadays, at lower temperature, the cross section 
can be the same for the s-waves or only smaller in 
case of p-waves annihilation channel [9]. The DM 
annihilation (DMA) signal is decreasing fast with 
the DM mass, at least as q ~ m~ 3 , but the p max 
is not well known and the significant DMA signal 
still can be obtained even at large masses of DM 
particles. Since the primordial clumps are much 
denser than the bulk component, most of the an- 
nihilation signal will come from the core of most 
abundant smallest clumps. The contribution from 
clumps is usually expressed as a 'boost factor' 
b ~ ^'""T and b » 1, probably except galac- 
tic center for the cusped bulk propfile. Taking a 
clump with M c i = 1O _8 M and 100 AU size with 
the average density of 100 GeV/cm 3 , the total yield 
for GeV charged particles will be <~ 10 23 s _1 for 
m x = 100 GeV. Assuming the isothermal distri- 
bution of clumps in the galactic halo of 20 kpc di- 
ameter and normalizing at n c i(8.5kpc) = 10pc~ 3 , 
the total luminosity from the DM annihilation in 
GeV range in 100 years will be <~ 10 45 particles, 
to be compared with the SNR explosion deliver- 
ing ~ 10 51 particles in the galactic disk. The DM 
clumps is a compact and constant source of CR and 



( Pcr ) ~ 10" 1L 

density distribution. 



3 producing a gradient in CR 



MHD turbulences initiated by DM an- 
nihilation 

The streaming of charged DMA products from the 
clump with the drift velocity above Alfven speed 
can increase the level of local MHD turbulences. 
The amplification of MHD waves parallel to mag- 
netic field lines can be calculated as [2]: 

G{r,k) « ^ffdpd»vp 2 (l - p 2 )5(p\p\ - 
ir) x {if + —if) 

he ' v a/i v op ' 

,where p is the cosine of scattering angle and 

%L ~ A^e^w % is the anisotropic term of the CR 
density distribution which can be obtained from 
integration of diffusion equation (2). The trans- 
verse waves are averaged out along propagation 
path and scattering is not efficient, that is, only tur- 
bulences along local field lines will be effectively 
amplified by streaming. The growth is reduced in 
dense molecular clouds with tig ~ 10 — 10 5 cto~ 3 
due to ions-neutral friction which dissipates en- 
ergy as S H ~ \(ocoiv a )n H , where (a co iv) ~ 
10 _9 cm 3 s _1 is the collisional cross section. The 
collision of opposite waves leads to a nonlinear 
damping proportional to the level of turbulences 
Sni « ^Wrf = S° nl W, where v t is the gas ther- 
mal velocity [3]. In the absence of external source 
of MHD waves the simplified solution of equation 
(1) and (2) neglecting energy dependencies, can be 
calculated for the spatial component analytically: 



W(r) 



exp(—g/r — s H r) 



r 2 (C + s n iexp(-g/r - s H r)/g) 



(3) 



,where g ~ {cv)Y tot ^f- is related to the total 

clump liminosity, s H = Sn/va, s n i = S° t /v a « 
sh are the dampings, and Co is a normalization 
factor. In the steady state it will result in fast in- 
crease of W(r) near the clump core followed by 
an exponential decrease due to damping, see Fig- 
ure 1 . The CR density distribution, obtained from 
solution of equation (2) with the diffusion defined 
by (3), will follow the W(r) dependency. The 
source function for spectral density W(k) is cutoff 
at k ~ ^- but can be extended to smaller wave 
numbers by collisions, reproducing Kolmogorov 
spectrum. More detailed consideration would in- 
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elude energy losses, fragmentation and interplay 
between electrons and protons streaming which is 
out of the scope of present paper. For the sta- 
ble growth the velocity of the clump proper mo- 
tion and the convection speed should be below 
local v a . In this case the anisotropic streaming 
of DMA products can create a region around the 
clump with small diffusion D c \ < D ext , aligned 
with the local B field and with size defined by the 
DM density in the clump, annihilation yields and 
dumping rates. The growth will be strongly sup- 
pressed inside dense gas clouds although particles 
produced in the cusp of the clump or on the bor- 
der of the cloud still can be confined. The energy 
losses and fragmentaion in the confinement region 
r c f can modify spectra of annihilation products. 
The time scale for the energy losses t; oss (E) has 
to be compared with the confinemnt time T c f ~ 
rl f /6D c i(E). The DMA antiprotons(protons) will 
lose energy in the gas clouds by ionization below 
100 MeV and by nuclear interactions at higher en- 
ergies. For positrons(electrons) the synchrotron 
and inverse Compton are the dominant losses at 
high and the bremsstrahlung at low energies. The 
modification of the energy spectra from the clump 
due to losses is shown in Figure 2 for different 
diffusion coefficients. Thus the contribution of 
charged DMA products to the locally observed CR 
fluxes will depend upon local environments and lo- 
cation of nearest clumps. 

Gamma rays from the DM clump will have two 
components, see Figure 3 where gamma rays have 
been calculated with the GALPROP code [8] mod- 
ified to include DMA and small scale variations in 
propagation parameters. First, the direct gamma 
rays from DMA \X ~ > 7 + X will produce a 
bump in the spectrum at E ~ 0.1 m x indepen- 
dent on the confinement [9] . Second, the elec- 
trons(positrons) will contribute to the lower en- 
ergy gamma rays via bremsstrahlung and inverse 
Compton, and protons(antiprotons) to the ir° peak 
but with somewhat harder spectrum as compare 
with the gamma rays from CR. The synchrotron 
radiation of electrons and positrons will also con- 
tribute to the radio waves. This can be an important 
source of radio emission in the galactic halo where 
the MHD enhancement can result in 5B/B > 1. 
The secondary contributions will strongly depend 
on the MHD enhancement and environments: gas 




r, AU 



Figure 1: The normalized W(r) and the f(r) 
density of GeV protons in the DM clump 

(Pmax ~ 10 3 GeV I cm 3 , m x = 100GeV,v a ~ 
lQ 6 cm/s,B ~ 1/zG) 

clouds, low density regions in the galactic disk 
or the galactic halo. The spectral features of the 
gamma DMA signal can be distinguished from the 
diffusive gamma rays produced by ordinary CR 
[10]. 

The confinement regions will also affect propaga- 
tion of the galactic CR. The low energy external 
CR will not penetrate deep inside such a clump, 
their contribution will be suppressed while the con- 
tribution of DMA particles produced in the clump 
will be enhanced. 

Conclusion 

The streaming of charged particles from DM an- 
nihilation in the cuspy DM clump can locally en- 
hance level of MHD turbulences reducing the lo- 
cal diffusion coefficient by orders of magnitude. 
The size of confinement region depends on the 
luminosity of the clump and damping of turbu- 
lences in dense gas clouds. The confinement re- 
gions will contribute to the small scale variations 
of propagation parameters and therefore CR den- 
sity. The gamma rays from the DM annihilation in 
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Figure 2: Modified by energy looses spectra of 
antiprotons and positrons from a clump(r c / = 
10- V) at different D cl /D ext . 



Figure 3: Diffusive spectrum of gamma rays from 
CR in conventional model(CM) [8] and the contri- 
bution from DMA (m x = 70GeV). 



the clump can be observed as a point like source 
with a particular spectrum with the bump from the 
direct gamma production in DM annihilation and 
secondary gamma rays produced by the confined 
charged particles in the clump. The annihilation in 
DM clumps will also produce synchrotron radia- 
tion even far away from the galactic disk. The con- 
tribution of charged DMA products to the observed 
CR fluxes will depend on the local environment 
and can be modified by energy losses. This can 
change relations between gamma rays and charged 
components from DMA. 
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